Environment-governed dynamics in driven quantum systems 
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We show that the dynamics of a driven quantum system weakly coupled to the environment can 
exhibit two distinct regimes. While the relaxation basis is usually determined by the system+drive 
Hamiltonian (system-governed dynamics), we find that under certain conditions it is determined 
by specific features of the environment, such as, the form of the coupling operator (environment- 
governed dynamics). We provide an effective coupling parameter describing the transition between 
the two regimes and discuss how to observe the transition in a superconducting charge pump. 



Introduction. Understanding how quantum systems in- 
teract with the environment [T] is of paramount impor- 
tance in quantum information science. While unveiling 
how the classical world emerges from the quantum one 
[2], it can also lead to a better protection against deco- 
herence effects on the way towards the realization of a 
quantum computer [3J. 

A standard approach to the dynamics of open quan- 
tum systems boils the problem down to the measurement 
of decohcrcncc rates, distinguishing between coherence 
loss, or dephasing, and relaxation. While this approach 
has successfully described a variety of quantum systems, 
it only offers a limited insight into the dynamics of de- 
coherence. A promising line of work developed in the 
last decade exploits the possibility of coupling the sys- 
tem to an engineered reservoir [4 12 . The chief advan- 
tage of engineered reservoirs is their tunability, allow- 
ing one to change at will environment features such as 
the spectral density, the temperature and the system- 
environment coupling strength. 

As new and more accurate ways are found of har- 
nessing the dynamic evolution of quantum systems, it 
becomes increasingly important to understand how the 
interaction with the environment is affected by a time- 
dependent modulation of the system parameters. Indeed, 
the study of dissipation in driven quantum systems is 
a long-established topic [13] that keeps finding new ap- 
plications, from quantum pumping |14Hlb| to quantum 
computation [TTlfTB] . 

In this Letter, we consider a periodically driven quan- 
tum system in the presence of a weakly coupled envi- 
ronment. We show that under certain conditions deco- 
herence takes place in a preferred basis determined by 
specific features of the environment, such as, the type 
of noise, rather than of the system. We label this un- 
usual regime as environment-governed dynamics (EGD), 
as opposed to the more familiar system-governed dynam- 
ics (SGD). We introduce an effective coupling parameter 
that presides over the transition between SGD and EGD. 
This parameter can be tuned by changing the properties 
of the drive. Our analysis is general and applies to opti- 



cal and solid-state systems alike. As a relevant example, 
we propose to observe the transition in a superconduct- 
ing charge pump [19, 20 . In this system the transition is 
controlled by an accessible experimental parameter and 
can be explored by measuring the pumped charge. 



Floquet- Born- Markov master equation. We consider 
a quantum system whose unitary evolution is governed 
by a periodic Hamiltonian H, so that H(t) = H(t + t), 
where r is the period. According to Floquet theorem, the 
Schrodinger equation admits solutions (Floquet states) 
of the form |* a (t)) = e~ zeat / h \ip a (t)} , where the Flo- 
quet mode \ip a (t)} satisfies \(p a (t + t)) — \ip a (t)) and 
e a is its corresponding quasienergy. Quasienergies and 
their associated modes are defined up to the transla- 
tion e a — > e a + Ml, where Q, = 2tt/t. As such, all 
quasienergies can be mapped into the first Brillouin zone 
[— gftfl, iM2]. The Hamiltonian describing the system 
and its environment is given by H to t = H(t) + HE + HsE, 
where He describes the environmental degrees of free- 
dom and Hse is the interaction term, that we assume of 
the form H$e = g S ® E, where g is the adimensional 
system-environement coupling constant and S and E are 
operators acting in the Hilbert space of the system and 
the environment, respectively. 

The general procedure for deriving the master equa- 
tion (ME) in the Floquet basis and in the Born-Markov 
approximation is outlined in Refs. 13 and 21. The 
resulting superoperator describing the time evolution 
of the density matrix is expressed a series of time- 
independent coefficients multiplied by phase factors of 
the form e lAa - l3 - k ~ iA ~<- s < k ' , with A. a ,p,k = + kCl. 

Starting from the general expression, we then perform 
a partial secular approximation (PSA). This amounts to 
neglecting all terms with k =/= k' while keeping those with 
k = k 1 and e a ^ ep. A residual time dependence due to 
terms oscillating like e^ e ° _e ' 3 ^ disappears when passing 
from the basis of Floquet states to that of Floquet modes. 
Our result, written in the basis of Floquet modes and in 
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the Schrodinger picture, reads: 
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We have introduced the following quantities: 



S(w) = J(w)no,(w) + 0(-w)J(-w) [1 + »Hh(-w)] 
A a/3fc = e Q - + fcfi , 

^e- lfcnt ((^ Q |B| W ) , 
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where 9(u>) is the Heaviside function, J(ui) the spectral 
density of the bath and nth(w) is the Bose-Einstein dis- 
tribution. 

A few remarks are in order. First, cannot be cast 
into a Pauli ME [25], i.e., the equations for the diagonal 
and off-diagonal terms are still coupled. This is due to 
the fact that we performed a PSA instead of a full secular 
approximation. The contribution of these terms becomes 
important close to degeneracies in the Floquet spectrum, 
where a full secular approximation is expected to fail [23 
121] . The PSA itself is justified provided the drive period 
is much faster than the decoherence time. This may not 
be true in the adiabatic limit, where other terms should 
instead be retained [23] , 

Equation ([TJ is formally akin to that of an undriven 
system, with quasienergies playing the role of ordinary 
eigenenergies and effective rates given by (§. The pres- 
ence of a drive manifests itself in the sum over k in Q , 
allowing the system to exchange energy with the environ- 
ment in any of the amounts A a g f.. The magnitude of 
this exchange is determined by the noise matrix elements 
X a ,p,k, telling whether the Floquet modes possess those 
energies and to what extent the coupling operator S al- 
lows the energy transfer to take place. Note that while 
each individual term in the sums ^ satisfies the detailed 
balance, the overall rates in general do not. 

Environment-governed dynamics. Equation (JIJ sug- 
gests that the dynamics of the coherences is determined 
by a competition between two terms. Let us now focus 
our attention on two states, a and j3. For a f3, the 
first term in ([lj is the Floquet energy gap £ — e a — eg, 
stemming from the nondissipative dynamics of the driven 



quantum system. The second term describes the effect 
of dissipation. For sufficiently small values of £, that 
is, close to a degeneracy in the Floquet spectrum, the 
dissipative term can dominate in (JT|). As a result, the 
dynamics of states a, (3 is strongly affected by the envi- 
ronment even if the system and the environment are only 
weakly coupled, i.e., g <C 1. This is due to the presence 
of a nearly resonant driving field introducing an energy 
scale, the Floquet gap, that can be much smaller than 
those of the undriven system. 

Let us now estimate the magnitude of the rates r in 
([3]). We consider an Ohmic environment, whose spectral 
density is given by J(uj) = aw/ c (w/w c ), where a is a 
dimensionless constant, f c a cutoff function and uj c the 
cutoff frequency. In the limit £1 <C uj c our results are 
independent of lo c and the explicit form of f c . As the 
sum in (J3J) runs over k, we may expect the dominant 
contributions to come from high-frequency modes (large 
k) such that X a gk does not vanish and kfl < u) c . These 
contributions are of order g 2 (£ + kCl) rs kg 2 Q,. 

We define an effective coupling parameter as 



x = g 2 n/£ 



(•5) 



If x "C 1 we are in the SGD regime. This corresponds 
to a decoherence time much longer than the period of 
Rabi oscillations between Floquet modes. By contrast, 
when x 1 we enter the EGD regime. This regime is 
realized close to a (quasi)degeneracy of the Floquet spec- 
trum, where Rabi oscillations between Floquet modes are 
slower than the decoherence time. We remark that the 
coupling parameter x is controlled by the Floquet gap £ , 
which can be tuned by changing the drive parameters. 
To appreciate the key role played by the drive, we ob- 
serve that an undriven system with an arbitrarily small 
energy gap cannot exhibit EGD. This is due to the fact 
that for any Ohmic and super-Ohmic environment the 
transition rates vanish for a vanishing energy gap. 

SGD to EGD transition. We will now provide an ex- 
plicit example of SGD to EGD transition in a driven 
two-level system. Using a pseudo-spin formalism, we 
write the system Hamiltonian as H(t) = a ■ B(t), where 
^ = Wx, Cyi o~z} ar e the Pauli operators and B(t) an 
effective magnetic field. Specifically, we consider a drive 
that modulates B along the loop shown in Fig. Qa). This 
drive is a realization of Landau-Zencr-Stuckclberg inter- 
ference [2S] with geometric phases [27]. While most of 
our results are far more general and can be obtained even 
using a single monochromatic drive, the geometric prop- 
erties of this particular drive provide a convenient way of 
tailoring the Floquet spectrum. In the following, we will 
study the effects of varying the angle ip [see Fig. [TJa)], 
which naturally controls the solid angle spanned by B in 
pseudo-spin space during a drive period. 

In Fig. [IJb) we plot the quasienergy spectrum of 
the system versus ip. The plot is obtained by numeri- 
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FIG. 1. (Color online) A model driven two-level system, (a) 
Loop drawn by the effective magnetic field B in pseudo-spin 
space. The solid angle spanned by B is controlled by the angle 
p. (b) Quasienergy spectrum as a function of ip. A weakly 
avoided quasienergy crossing occurs at p c = —0.26. The drive 
parameters are the same as in Fig. [3]with the definitions given 
in the text. 



cally solving the Schrodinger equation for the evolution 
operator generated by H{t). The quasienergy gap E 
sharply decreases near tp = tp c , where a weakly avoided 
quasienergy crossing occurs. Close to cp c , we thus expect 
EGD to be attained. 

We characterize the transition from SGD to EGD by 
studying the quasistationary state approached by the 
driven system in the presence of dissipation. For the 
purpose of illustration, we consider a zero-temperature 
environment. Our results are obtained by enforcing 
Pa/3 = in (JTJ. This results in a set of algebraic equa- 
tions for p with coefficients X a ^^ and £. Their solu- 
tion gives the quasistationary density matrix p st , which 
is time-independent in the Floquet-mode basis. We write 
p st = 1/2(1 + n-S) where £ is the vector of Pauli oper- 
ators in the Floquet-mode basis. In this way, the resid- 
ual coherence between Floquet states is associated to the 
quantity n± - 



(c) 



In Fig. [2^a) we plot n z (full line) and nj_ (dashed line) 
versus (p in a neighborhood of ip c . We choose the noise 
operator S = a z and g 2 = 0.01. As ip approaches <p c from 
either side, we witness the transition from SGD to EGD, 
signalled by a revival of the coherence between Floquet 
states. 

We then fix tp to a value close to tp c and change x 
by changing g 2 , to which \ is proportional. The re- 
sults are presented in Fig. |2jb). For small values of g 2 , 
the coherence between Floquet states is completely lost 
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FIG. 2. (Color online) SGD to EGD transition. n z (full lines) 
and iii (dashed lines) versus <p for S = a z and g 2 = 0.01 (a), 
versus g 2 for S = a z and tp = —0.24 (b), and versus a for 
ip = -0.16 (black) and ip = -0.25 (red) and g 2 = 0.1 (c). 
EGD is attained is the highlighted regions of panels (a,b) and 
everywhere in panel (c). The drive parameters are the same 
as in Fig. [3] 



(nj_ = 0, SGD). For larger values of g 2 (but still in the 
weak-coupling regime), n± takes a finite value (EGD). 
Notice that as soon as either of the two limits is attained, 
the steady-state populations do not depend on the exact 
value of g 2 . 

A notable feature of the EGD regime is that the pre- 
ferred basis for relaxation is strongly affected by the type 
of noise. To demonstrate this, we introduce a family 
of coupling operators 5(a) = aa x + (1 — a)a z , where 
a G [0, 1] parametrizes the angle between the reference 
basis of H and the noise operator. In Fig. [2^c) we then 
plot n z (solid lines) and nj_ (dashed lines) versus a for 
two different values of ip. We set g 2 = 0.1, a value ensur- 
ing that the EGD limit is attained for both cases. Upon 
changing a, p s t undergoes significant changes. This is in 
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FIG. 3. Application to the Cooper-pair sluice, (a) Equivalent 
circuit of the sluice, (b) Time dependence of Jl, Jr and 
n g during a pumping cycle, (c) Pumped charge Q p versus ip 
for selected values of g 2 . The relevant drive parameters [see 
panel (b)] are: r=l ns, Ec=l K, An g = 0.2, J max = O.lEc, 
7—1 fl~ 3 J 



stark contrast to what happens in SGD, where relaxation 
always takes place in the same (Floquet) basis regardless 
of the coupling operator. Note, however, that even in 
SGD different coupling operators can result in different 
populations, as the detailed balance is not satisfied. 

As the quasistationary dynamics is determined by a 
set of algebraic equations, many of these features can 
be addressed analytically. In the SGD limit (v <c 1), 
we obtain the same results as predicted by a full secular 
approximation. We find pf* = pn,SGT>> where Ph,sgd is 
a constant whose explicit expression is found in literature 
[T51 US], and pf 2 = 0(x)- By contrast, in the EGD limit 
(X > 1) we obtain pf x = /9h,egd and pf 2 = Pi 2 ,egd + 
0(l/x). Explicit expressions for the constants Pi^sgd, 
Pii,egd and pi2,EGD are provided in the Supplemental 
Material. 

Altogether, the results of Fig.[2]show that in the EGD 
regime the environment is drastically influencing the dy- 
namics of the system. In principle, a steady-state mea- 
surement of «£c,y,z can thus disclose valuable information 
on both the type of noise and the strength of system- 
environment coupling. 

Observation in a superconducting charge pump. We 
now discuss how to observe the SGD to EGD transi- 
tion in a superconducting charge pump, the Cooper- 
pair sluice [TSJ [213 HI]- The sluice, shown schemati- 
cally in Fig. [3ja), consists of a single superconducting 
island, coupled to superconducting leads via two su- 
perconducting quantum interference devices (SQUIDs). 
The SQUIDs are operated as Josephson junctions of 



tunable energies JL,n(t). A gate electrode capacitively 
coupled to the island provides a third control param- 
eter by inducing a polarization charge n g (t) in units 
of 2e. The device is operated under a constant su- 
perconducting phase bias tps- In the charging regime 
Ec 3> Jl,b, (Ec is the charging energy of the island), 
the dynamics can be reduced to the two lowest-energy 
charge states |0) and |1) corresponding to zero and one 
excess Cooper pairs on the island, respectively. The 
system is then described by a pseudospin Hamiltonian 
of the form discussed above, with effective field compo- 
nents B x (t) = ^ J + (t) cos ^, B y (t) — J_ (t) sin ^ , and 
B z (t) = E c [1/2 - n g (t)}, where J±(i) = J L {t) ± J R (t). 
Pumping is achieved by steering the three parameters Jl , 
Jr and n g in a periodic fashion, as shown in Fig.|3|b). 

Starting from the given definitions, it can be shown 
that the pumping cycle of the sluice is a realization of 
the loop of Fig. [TJa), with ip — ips- The quasienergy 
gap of the sluice can thus be tuned by changing the su- 
perconducting phase bias while performing exactly the 
same pulse sequence. This makes the sluice an excellent 
candidate to verify our theoretical predictions. 

A further advantage in using the sluice is the possi- 
bility to read out the pumped charge Q p . The formal 
definition of Q p , its behavior in the adiabatic limit and 
its connection with the Berry phase are well understood 
[5TJ1 [51] . In realistic experimental conditions [29 ] 32] , the 
measured Q p is often averaged over a great number of cy- 
cles. As such, it reflects the quasistationary state reached 
by the system in the presence of dissipation. 

The main source of decoherence in the sluice is charge 
noise, due to fluctuations in the gate voltage [71Q3]. We 
describe it by putting S = a z , a = 1 and g = Cg/C^., 
where C g and Cs are the gate-to-island capacitance and 
total island capacitance, respectively. 

In Fig. [3^c) we plot Q p as a function of ip for differ- 
ent coupling strengths g 2 . For small values of g 2 , corre- 
sponding to x ^ 1 , Qp exhibits a dip around tp c . The dip 
appears in coincidence with the weakly avoided crossing 
in the quasienergy spectrum [Fig. [TJb)], and stems from 
the mixing of "adiabatic" Floquet states at the crossing 
[TSJ - As g 2 is increased, an expanding neighborhood 
of (p c undergoes the transition to EGD, producing an in- 
crease in Q p . Finally, for large enough values of g 2 the 
whole region of mixing is in the EGD regime, and the dip 
has disappeared. 

Conclusions. A driven quantum system interacting 
with the environment exhibits a richer scenario than an 
undriven one. This is due to the emergence of an energy 
scale, the Floquet gap, that can compete with decoher- 
ence rates in the vicinity of a quasienergy crossing. This 
energy can be tuned by choosing the drive parameters. 
We have identified two dynamical regimes and given an 
effective coupling parameter governing the transition be- 
tween the two. This transition manifests itself in the 
quasistationary density matrix, in particular, in the re- 
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vival of coherences between Floquet states. Finally, we 
have discussed how to observe the transition in a super- 
conducting charge pump. A closely related system that 
can be considered in a similar spirit is the driven Cooper- 
pair box [551155]. 

While the SGD regime has been intensively studied, 
the EGD regime is vastly unexplored. Due to its simplic- 
ity and applicability to a variety of different systems, the 
present approach may emerge as a useful tool to study 
system-environment interactions in open quantum sys- 
tems. 
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Quasistationary solution for the master equation 



Starting from Eq. ([I]) in the main text, we restrict ourselves to a two- level system and set p a p = 0. We arc not 
here interested in the details of the solution but only in its scaling with respect to the quantities g, £1 and £. 

We represent the quasistationary solution p st as a vector of components (pf[ , pf 2 , p\ \ , P22 ) ■ Using (JTJ , the definitions 
([2]) and (|3| and the symmetries of the rates pj, we find that p st satisfies the following matrix equation: 
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where m, . . . ,s are a set of dimensionless coefficients. We have factored out a factor g 2 £l in the Redfield tensor ^ as 
this is its expected scaling in the case we consider (Ohmic environment, <^ uj c ). 

For simplicity, let us assume that all coefficients in ^ are real (it turns out that this is often the case). The exact 
solution is then given by: 
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with 



(7) 



A = (q - r)[o{q + r) - 2ns]g 4 n 2 + 0Z 2 , 
C = 2n(q -r)(p- s)g A n 2 - (m - o) [{q 2 
D = (ms - op)[(q - r)g 2 il + i£]g 2 n . 



In the SGD regime x < 1, A m o£ 2 , C ■■ 



(o — m)£ 2 and D « i(ms 
o 



r 2 )g^ 2 + £ 2 )} , 

■ op)g 2 fl£. An approximate solution is 
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Pll, SGD , 

= 0(x) , 
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(8a) 
(8b) 



The coherence pf^ vanishes for infinitesimally weak coupling, indicating that the Floquet basis is the relaxation basis. 
Furthermore, (8a) does not depend on the coupling strength. These result are the same as obtained with a full secular 
approximation |13| 115] . 

In the EGD regime \ ~ 1> an approximate solution is: 



P\l 



P S 12 



o{q 



2ns 



2n(p — s) — (to — o)(q + r) 

ms — op 
2n(jp — s) — (to — o)(q + r) 
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(9a) 
(9b) 



The coherence pf 2 in general does not vanish; in the limit \ 3> 1, it also approaches a constant value. 



